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Abstract 

In recent years, a market for mortality derivatives began developing as a way to handle system- 
atic mortality risk, which is inherent in life insurance and annuity contracts. Systematic mortality 
risk is due to the uncertain development of future mortality intensities, or hazard rates. In this 
paper, we develop a theory for pricing pure endowments when hedging with a mortality forward 
is allowed. The hazard rate associated with the pure endowment and the reference hazard rate for 
the mortality forward are correlated and are modeled by diffusion processes. We price the pure 
endowment by assuming that the issuing company hedges its contract with the mortality forward 
and requires compensation for the unhedgeable part of the mortality risk in the form of a pre- 
specified instantaneous Sharpe ratio. The major result of this paper is that the value per contract 
solves a linear partial differential equation as the number of contracts approaches infinity. One can 
represent the limiting price as an expectation under an equivalent martingale measure. Another 
important result is that hedging with the mortality forward may raise or lower the price of this 



pure endowment comparing to its price without hedging, as determined in Bayraktar et al. 2009 
The market price of the reference mortality risk and the correlation between the two portfolios 
jointly determine the cost of hedging. We demonstrate our results using numerical examples. 

Keywords. Life annuities, longevity risk, g-forward, mortality-linked derivatives, instantaneous 
Sharpe ratio, incomplete market 



1 Introduction and Motivation 



A basic assumption in many actuarial texts is that mortality risk can be eliminated based on the 
law of large number. It is believed that the standard deviation per insurance policy vanishes as 
the number of policies sold becomes large enough. However, this assumption is valid only when the 
mortality intensity is deterministic, and a number of recent researchers argue that mortality intensity, 



or hazard rate, is stochastic; see, for example, Dowd et al. 2006 and the references therein. The 



uncertainty of hazard rates is significant enough that stochastic mortality risk has to be considered 
in the valuation of life insurance and annuity contracts and in pension fund management. A concrete 
example of stochastic mortality risk is longevity risk, namely, the risk that future lifetimes will be 
greater than expected. Longevity risk has attracted much attention in recent years, and many capital 
market instruments have been proposed to deal with this risk for annuity providers and pension funds; 
see 



Dowd et al. 2006 , Blake and Burrows 2001 , and Blake et al. 2006 for more details 



However, few researchers have focused on the effectiveness of hedging mortality risk with the 

In 



proposed mortality- linked derivatives; one notable exception is the work of Lin and Cox 2005 
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this paper, we investigate the apphcation of mortahty-hnked derivatives for hedging mortahty risk 
and offer suggestions for further mortahty-hnked innovation based on our analysis. To this end, we 
select a stochastic model to describe the mortality dynamics. Several stochastic mortality models 



2001 , 


Biffis 


2005 , 


Schrager 



2006 , and Dahl 2004 use continuous-time diffusion processes to model the hazard rate, as we do in 



this paper. Alternatively, Miltersen and Persson 2005 and Cairns et al. 2006a model the forward 
mortality. Milidonis et al. 2010 incorporate mortality state changes into the mortality dynamics with 
a discrete-time Markov regime-switching model. Also, see Cairns et al. 2006b for a detailed overview 



of various modeling frameworks. In this paper, we use the model proposed by Bayraktar et al. 2009 
to describe the dynamics of both hazard rates: A^, the one inherent in the insurance contract to be 
hedged, and A^, the one referenced by the mortality-linked derivative. 

Another issue is the choice of pricing paradigm. Different methods for pricing mortality risk have 
been proposed in recent literatures, and Bauer et al. 2010 extensively discusses them. Among these 
methods, Bayraktar et al. 2009 developed a dynamic pricing theory, which can be considered as a 
continuous version of the actuarial standard deviation premium principle. In our paper, we extend 
their pricing mechanism to a market that includes mortality-linked derivatives. We price a pure 
endowment assuming that the issuing company hedges its contract with a mortality forward in order 
to minimize the variance of the value of the hedging portfolio and then requires compensation for the 
unhedgeable part of the mortality risk in the form of a pre-specified instantaneous Sharpe ratio. 

The main purpose of this paper is to investigate the hedging of life insurance and annuity contracts 
with mortality-linked derivatives. To this end, we develop a partial differential equation (PDE) whose 
solution is the value of the hedged insurance contract. We compare the values of the hedged contract 
under different market prices of mortality risk. We also analyze how the correlation between Af and Xj 
affects the values of the hedged contract. The main contribution of our paper is to show that hedging 
can reduce the price of the insurance contract only under certain conditions on the correlation of the 
hazard rates and on the market price of mortality risk. As part of the procedure, we also show that 



the desired features of the pricing mechanism by Milevsky et al. 2005 and Bayraktar et al. 2009 
still hold in our extension. 

The remainder of this paper is organized as follows: In Section[2] we present our financial market, 
describe the pricing mechanism of the pure endowment in a market with mortality-linked derivatives, 
and derive a non-linear PDE whose solution is the value of the hedged pure endowment. In Section 
[sj we analyze the value p(") of n pure endowments on conditionally independent and identically 
distributed lives, with the emphasis on how the correlation of the hazard rates and the market price 
of mortality risk affect the price of the hedged pure endowments. We then present the PDE that gives 
the limiting value of ^-P*-"' as n goes to infinity in Section |4| We show that this limiting value solves 
a linear PDE and represent this value as an expectation with respect to an equivalent martingale 
measure. In Section [5| we demonstrate our results with numerical examples, discuss whether and 
when the hedging with mortality-risk derivatives reduces the price of pure endowments, and provide 
suggestions on the application of mortality-linked derivatives for insurance companies. We describe a 
numerical scheme to compute the value of a pure endowment in Section [7] Section [6] concludes the 
paper. 
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2 Incomplete Market of Financial and Mortality Derivatives 



In this section, we describe the pure endowment contract and the financial market in which the issuer 
of the contract invests to hedge the risk. In the financial market, there are three products: a money 
market fund, a bond, and a mortality derivative. We obtain the optimal strategy to hedge the risk of 
the contract with bonds and mortality derivatives in order to minimize the variance of the value of 
the investment portfolio. We, then, price the pure endowment using the instantaneous Sharpe ratio. 

2.1 Mortality Model and Financial Market 

First, we set up the model for the dynamics of hazard rates-either the hazard rate for the pure 
endowment or the one for the mortality derivative. We assume that a hazard rate At follows a diffusion 
process with some positive lower bound A. Thus, we require that as Aj goes to A, the drift of At is 
positive and the volatility of At approaches 0. Biologically, the lower bound A represents the remaining 
hazard rate after all accidental or preventable causes of death have been removed. Mathematically, 
the need for such a lower bound appears later in this paper. 

Specifically, we use the following diffusion model for a hazard rate: 

dA( = a(At, t) (At -X)dt + b{t) (At - A) dWt, (2.1) 

in which W is a standard Brownian motion on a filtered probability space (n,F, (Ft)t>0)P)- We 
require that the volatility b{t) is a continuous function of t and is bounded from below by a positive 
constant k in [0,T]. We also assume that a(Af,t) is Holder continuous with respect to A and t, and 
that a(At, t) > when < At — A < e for some e > 0. 

In this paper, we consider two different but correlated hazard rates. One is the hazard rate of the 
insured population; namely, the hazard rate of the people who purchase the pure endowments. For 
simplicity, when we consider a portfolio of n pure endowment contracts in this paper, we assume that 
all individuals are of the same age and are subject to the same hazard rate. We denote as Af the 
hazard rate of insured population, and the dynamics of \f is given by 

dAf = a^(Ar, t) (Af - A^) di + h^{t) (Af - A^) dWf . (2.2) 

We also consider a second hazard rate on which the mortality derivatives are based, namely, the 
hazard rate of an indexed population. We denote this hazard rate as A^ , whose dynamics is given by 

dA,^ = a\\i, t) (\i - \') dt + h\t) (\{ - \') AWi (2.3) 

The uncertainties of the two hazard rates are correlated such that dW/ dW/" = pdt with p e [—1,1]. 

Suppose, at time t = 0, an insurer issues a pure endowment to an individual that pays $1 at time 
T if the individual is alive at that time. To price this contract, we will create a portfolio composed of 
the obligation to pay this pure endowment and investment in the financial market. 

In the financial market, the dynamics of the short rate rt is given by 

drt = /x(rt, t) dt + a(rt, t) dW^ (2.4) 

in which fj, and cr > are deterministic functions of the short rate and time, and is a standard 
Brownian motion adapted to (n,F, (Ft)t>o,P)- We assume that is independent of and W\ 
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and that /i and a are such that rf > almost surely for all i > and such that (2.4) has a unique 
solution. 

Both the T-bond and the mortality derivative are priced based on the principle of no-arbitrage. 
Thus, for the short rate r, there exists a market price of risk that is adapted to the filtration 
generated by W^; and for the hazard rate Aj , there exists a market price of risk that is adapted 
to the filtration generated by . We, therefore, write cither ql or q^'{rt,t) for the market price of 
the short rate risk at time t; similarly, we write either q^ or {^(t^) ^'^^ market price of the 
hazard rate risk at time t. 

Define an equivalent martingale measure Q whose Radon-Nikodym derivative with respect to P is 
given by 



dF 



exp ■ 



q^{r,,s)dW: + q'\xi,s)dW^]-^J^ \iq^ (r,, s)f + [q'' [Xl s)) 



In the Q-space, the dynamics of the hazard rates and the short rate are given by 



dAf 



a^.Q(Ar, A^, t) (Af - A^) dt + bP{t) (Af - A^) 
dA^ = a'^Q{\i, t) (\l ~ A^) dt + b'{t) (\i ~ \') dW^/''3 
drt = ^Ji'^{rut) dt + G{rt,t) dWl^^ , 



dW, 



P.Q 



ds 



(2.5) 



(2.6) 



in which 



and 



P,Q 



^wr+pj;;q^'{\i,s)ds, 



i,Q _ 

^ — 



Wi + r,q^ {\i,s)ds, 
Wl + jlq^{r,,s)ds, 



Af , t) = aP{\f,t) - pq^' {\l t) bP{t), 
a'Q{\lt) = a\\lt)-q^\\lt)b\t), 
^li'iirut) ^ p{rut) - q-{rut)a{rut). 

The time-t price of the T-bond is given by 



F{r,t;T) = 



- ds 



rt 



and the bond price F solves the following partial differential equation (PDE), Bjork 2004 : 



(2.7) 



(2. 



(2.9) 



Ft + Ai^(r, t) Fr + ^a^{r, t) F„ - rF = 0, 
F{r,T-T) - 1. 

From this PDE, we obtain the dynamics of F for t < s < T: 

dF(r„ s) = [r,F(r„ s) + g'-(r„ s) a(r„ s) F,(r„ s)] ds + a(r„ s) F,(r„ s) dWJ 
F{rt,t)^F{r,t). 



(2.10) 



(2.11) 



Without loss of generality, we specify the mortality derivative as a q- forward. Define the cumulative 
hazard rate process by Aj = Af ds for < i < T. Then, the time-t value of a g-forward with delivery 
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time T is given by 



rt 



r, A 



(2.12) 



in which K = EQ 



Fo 



is the dehvery price. As for the T-bond, we have the fohowing PDE 
for the mortahty derivative: 

St + + 5^x2 Srr + a^''^ ' (a^ - A^) 5-^/ + ^ {h'f (X^ - A^)' 5a/a^ + - 

S{r,\\A',T;T) = e-^' - K. 



rS = 0, 



From the PDE above, we obtain the dynamics oi S ioi t < s < T: 



dSs 



'sS + ql aSr + ql' ■ (a^ - A^) Sx' ds + a Sr dW,' + b' ■ (a^ - A^) Sxi dW^ , 



(2.13) 



(2.14) 



Since we fixed the maturity T for both the bond and the g-forward in this paper, we drop the notation 
T when appropriate. 



2.2 Pricing the Pure Endowment via the Instantaneous Sharpe Ratio 
2.2.1 Recipe for Valuation 

Even with the mortahty derivatives, the market for insurance is incomplete when p ^ ±1 due to the 
fact that the mortality of the insured population and the mortality of the indexed population are not 



perfectly correlated. This mismatch is called basis risk; see Coughlan et al. 2007 for more details 



Therefore, there is no unique method of pricing for insurance contracts, and to value contracts in this 



market, one has to choose a pricing mechanism. For example, Bayraktar and Ludkovski 2009 used 



indifference pricing and Dahl and MoUer 2006 consider the set of equivalent martingale measures 



when pricing the unhedgeable mortality risk. We use the instantaneous Sharpe ratio proposed by 



Milevsky et al. 2005 and Bayraktar et al. 2009 to price the risk due to the desirable properties 



of the resulting price. We will show these properties in Section [3] Moreover, as the number of 
contracts approaches infinity, the limiting price per contract can be represented as an expectation with 
respect to an equivalent martingale measure. In a market without mortality derivatives, this pricing 
methodology has been proved useful for pricing pure endowments ( Milevsky et al. 2005| ) , life insurance 
(Young [2008 ), life annuities (Bayraktar et al. 2009| ), and financial derivatives (Bayraktar and Young 



2007b| ). In this paper, we extend this pricing mechanism to incorporate mortality derivatives in the 



financial market. Our method for pricing in an incomplete market with mortality derivatives is as 
follows: 

1. First, we set up a portfolio composed of two parts: (1) the obligation to underwrite the pure 
endowment, and (2) a self-financing sub-portfolio of T-bonds, g-forwards maturing at T, and 
money market funds to partially hedge the pure endowment contract. 

2. Second, we find the optimal investments in bonds and mortality derivatives to minimize the local 

In 



variance of the portfolio. This method is called local risk minimization by Schweizer 2001 



case of a complete market, the minimized local volatility is zero. However, the incompleteness 
of the insurance market leads to residual risk, as measured by the local variance. 
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Third, we assume that the insurance provider requires compensations for the unhedgeable risk. 
The price of the contingent claim is set to make the instantaneous Sharpe ratio of the total 
portfolio equal to a pre-specified value. This is equivalent to setting the price of the contingent 
claim such that the drift of the portfolio equals the short rate times the portfolio value plus the 
pre-specified Sharpe ratio times the local standard deviation of the portfolio. Thus, our pricing 



method is a type of local standard deviation premium principle. Young 2004 



2.2.2 Hedging and Pricing a Single Pure Endowment 

Denote by P{r, , ,t;T) the time-i value of a pure endowment that pays $1 at maturity T if the 
individual is alive at that time. Here, we explicitly recognize that the price of the pure endowment 
depends on the short rate r, the hazard rate of the insured individual, and the hazard rate A^ of 
the indexed population. Since the maturity T is fixed, we simplify the notation to P{r, X^ ,X^ ,t). (By 
writing P to represent the value of the pure endowment, we assume that the individual is alive. If the 
individual dies before T, the value of the pure endowment jumps to $0.) 



Suppose the insurer creates a portfolio 11 as described in Step (i) in Section 2.2.1 This portfolio 



is consist of two parts: (1) the obligation to underwrite the pure endowment with value — P, and 
(2) a self-financing sub-portfolio Vt of T-bonds, g-forwards, and money market funds to hedge the 
risk of the pure endowment. Thus, Ilf = —P{rt,Xj,Xf,t) + Vt- Let ttJ" equal the number T-bonds 
and TT^ the number of g- forwards in the self- financing sub-portfolio at time t with the rest, namely, 
Vt — 7r[F(rt,t) — TT^' S{rt, X(,t), in money market funds. 

By Ito's lemma, the dynamics of the value of the pure endowment P{r, A^, X^ ,t) in the physical 
probability space is given by 



dP(r, A^ A^, t) =Ptdt + Pr drt + P^i dX^ + P^p dAf + -P^ d[r, r]t 

+ ^P^^^P d[X', X']t + P^i^p d[A^ X% + ^P^P^P d[A^, X% - P dNt 



Pt+f^Pr 
^Cj'Pr. 



K - A' 



2 



A^ ) P,P 



dt 



p. 



(Af 

2 \ {h^f{xl-X^)"p,.,. + \{b^f{x 

ph'h^ {K-h') (Ar-A^)Pv.pdt- 

a Pr dWl + ■ {x{ - A^) P^i dWl + 6^ • f Af - A^ ) P^p dW, 



(2.15) 



dt 



in which [•, ■]t represents the quadratic variation at time t, and Nt is a time-inhomogeneous Poisson 
process with intensity Af that indicates when the individual dies. Recall that the value of P jumps 
to $0 when the individual dies; thus, we have the term —PdNt to account for this drop. 
Since the sub-portfolio Vt is self-financing, its dynamics are given by 



dVt = < dP(rt, t) + n^'dSt {rt, Af t) + n \Vt - < P(rt, t) - t:^' St {n, Af t) 



dt 



< g[ aFr + Trf a Sr + t:^' q^' b' ■ (a^ - A^) S^i + n Vt 

■( 



dt 



(2.16) 



TTj a Fr 



a dWl + i^^'h' ■ (X\ - A 5,. dWl, 



in which the second equality follows from equations (2.11 ) and (2.14), and we suppress the dependence 
of the functions on the underlying variables. 
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It follows from equations (2.15) and (2.161 that the value of the portfolio Ilt+h at time t + h ioi 
h > 0, given Ilf = 11, is 



n 



t+h 



n 



t+h 



t+h 



dV, 



U-J DP{rs,Xf,Xi,s)ds 

h 

irlqlaFr + TT^'qlaSr + n^' q^' b' ■ (a^ - A^j S^i 

rt+h 



ds 



(2.17) 



crFr + TT^ aSr - crPr dW, 



- J^^ 6^ • (Af - A^) F;,P diyf + ^ ^ P (dN, - Xf ds) + J^^ rsH^ds, 

in which S) is the operator defined on the set of appropriately differentiable functions on x (A^, oo) x 
(A^,oo) X [0,r] by 

2)w = - (r + A^) v + vt+ fivr + ■ (^X' - X/^ vxi + ■ (a^ - A^^ vxp + \^<J^Vrr 

+ 2 {b^f - A^)' «A^A^ +Pb'b^ (a^ - A^) (a^ - A^) v,r,P + \ (b^f (a^ - A^)\,.,P. 

(2.18) 

Note that the compensated counting process Nt — Jp Af ds is a (local) martingale . 

When we consider the single life case, the value of Pt becomes zero immediately after the individ- 
ual's death, so the value of the portfolio increase by P. If we consider the price P^"^ of n conditionally 
independent and identically distributed lives as in Section 2.2.3 the intensity of the counting pro- 
cess Nt is nXf at time t. As one of the n individual dies, the value of the portfolio increases by 
p(") _ p("-i), \Ye will consider P*^"^ later and continue with the single-life case for now. 
The second step as stated in Section 



2.2.1 



is to choose ttJ" and tt^ to minimize the local variance 



of the portfolio. To this end, we calculate the conditional expectation and variance of Tit+h at time t 
given lit = n. First, we define a stochastic process Yh for ft, > by 



pt+h fZ-tfi 

Yh = U- DP(rs,Af,Af,s)ds+ / r,n,ds 

Jt Jt 

■h 

< q: aFr + TT^' ql a Sr + tt^' q^' ■ ( Af - AM 5v 



rt+h 



(2.19) 



ds. 



Thus, E(nt+^|Ft) = E'"''^^'^^'*(Yh,), in which E'''^^^^^'* denotes the conditional expectation given 
rt ^r, X't ^ X', and Af = A^. From ([2J7]) and ( pl9| we have 

nt+,, = n + (nlaFr + TT^'aSr - (tP.) dW^J + b' ■ (^Xi - A^) (n^' S^r - Pa^) dW^ 

rt+h 

PxpdW,^ + J P{dN,~Xfds). 



Af-A^ 



(2.20) 
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It follows that 

Var[nt+,,|Ft] =E 



E 



t+h 



aFr + TTg aSr — (jPr I ds 



E 



2E 



y*^%6^6^ (Af - A^) (Af - A^) P,. [n^'s.s - P,s) ds + K j'^'^ X^'P'ds + o{h), 



(2.21) 

in which all the expectations are conditional on the information available at time t. Thus, the optimal 
investments in the g-forward and T-bond to minimize the local variance are given by, respectively, 



1 
1 



Px'+P 



Pr- U- 



b^-(Ar-A^ 



P 



Sr 



(2.22) 



(2.23) 



Equations (2.22) and (2.23) show that in the self-financing sub-portfolio, the g-forward is used to 



hedge the mortality risk in the pure endowment, and T-bonds are used to hedge the interest risk of 
the portfolio. Under this investment strategy, the drift and local variance of the portfolio become 



lim \ [E(nt+JFO - n] = -S^P + rn. 



and 



with 



1 



lini -Var [nt+,,|Ft] = (1 - p") {b' 



A^-A' 



P 



A J' 



\PP\ 



(2.24) 
(2.25) 



S)«P = - (r + A^) P + P, + M«P, + a^'« f A^ - A^ P,. + a^-« f A^ - A^) P,p + ^a^P, 



+ I (h^f (a^ - P,.,. +pb'b^ (a^ - V) (a^ - A^) P,:,. + \ {b^f (a^ - A^)' P,P,. 

(2.26) 



Remark 2.1 When p = ±1, the q-forward and the pure endowment bear identical uncertainty risk in 
the hazard rates. In this case, the mortality risk in the pure endowment can be completely hedged with 
the q-forward, and the minimum local variance of the portfolio only comes from the random occurrence 
of death, namely. 



lim —Var 

h^O h 



t+h 



Ft 



A^p2. 



(2.27) 



Remark 2.2 As we will show in Property 3.8 in Section^ P\i = when p = Q, and the corresponding 
optimal investment in the q-forward is tt^ =0. Intuitively, the q-forward is not used to hedge the 
mortality risk in the pure endowment when the two underlying hazard rates are not correlated. 



Next, we price the pure endowment via the instantaneous Sharpe ratio as stated in Step (iii) in 
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Section 2.2.1 The minimized local variance of the portfolio in (2.25) is positive; therefore, the insurer 



is not able to hedge all the risk underlying the pure endowment. The insurer requires an excess return 
on this unhedgeable risk so that the instantaneous Sharpe ratio of the portfolio equals a pre-specified 
value a. We could allow a to be a function of say r, A''^, A^, and t to parallel the market price of the 
risk process |. However, for simplicity we choose a to be a constant. (Further discussion of 



the instantaneous Sharpe ratio is available in 



Milevsky et al. 



2006 



.) We assume that < a < 



as we will see, some of the properties of P rely on this upper bound for a. 

To achieve a Sharpe ratio of a and thereby to determine the value P of the pure endowment, we 
set the drift of the portfolio equal to short rate times the portfolio value plus a times the minimized 



local standard deviation of the portfolio. Thus, we get the following equation for P from (2.24) and 



(2.25) 



(2.28) 



If the individual is still alive at time T, then the policy is worth exactly $1 at that time, that is, 
P (r, A^, A^, T) = 1. Thus, P = P{r, \\ A^, t) solves the following non-hnear PDE on M+ x (A^, oo) x 
(A^,oo) X [0,r]: 



Pt + tiQPr + a^'Q ■ (a^ - A^) + a""'^ ■ (a 

+ la^P^r + \ {b'Y (a^ - A^)' P,.,. + i (b^Y ( 
+pb'bP (a^ - A^) (a^ - A^) P^i^P - (r + A^) P 

= (1 - p') {bn^ (a^ - A^) ' + \PP^ 

[P(r,A^A^,^) = l. 



P 



(2.29) 



We can simplify the solution to (2.29) because the uncertainty in the short rate is uncorrelated 



with the uncertainty in mortality rates. Indeed, note that P{r, X^ , X^ ,t) — F{r,t) tp (^X^ , X^ ,t) , in 



which F is the price of the T-bond and solves (2.10), and V' solves the following non-linear PDE: 

'^t + a^-Q • (a^ - A^) V.,. + a^-Q ■ (a^ - A^) + ^ {b'f (a^ - A^)%,.,. 
+pbV [X' X') (a^ - A^) + 1 {b^f (a^ - A^)%,P,P - A^V 



{1 - p^) {bPf (^XP - X"") yjl^+XP^2 
[V'(A^A^,T) = 1. 



(2.30) 



The existence of a solution to (2.30) follows from standard techniques; see, for example, Chapter 36 
in Walter 1970 . Uniqueness of the solution follows from the comparison result in Section [s] of this 



paper. 



2.2.3 Hedging and Pricing a Portfolio of Pure Endowments 

In this section, we develop the PDE for the price P^"^ of n pure endowment contracts. We assume 



that all the individuals are of the same age and are subject to the same hazard rate given in (2.2 1 



We further assume that, given the hazard rate, occurrences of death are independent. As discussed 



in the paragraph following equation (2.18), when an individual dies, the portfolio value 11 increases 
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by p(") — F^" ^\ By paralleling the derivation of (2.29), one gets the following PDE for 



>f"' + fiQp^''^ + a''Q ■ (a^ - A^) P^J^ + a^^Q ■ (a^ - A^) p|';) 

(a^ - A^) (a^ - A^) - rP(") - nA^ • (p(") - P("-i)) 

= -0-^(1 - p2) (feP)2 I^^P _ ^P^' {Px^^y + (P(") - P(«-i))', 
P("H7-,A-^,A°,T) = n, 



(2.31) 



with initial value P(°) = 0, and P'^^ = P, as given by ( |2.29 l 
As in Section 



2.2.2 



P(")(r, A^A^,i) = F{r,t)il)''-^\\^ ,XP ,t), in which P solves ([2lo| and V'^"^ 



solves the following PDE 



+ a^-Q • (a^ - A^) ^('J) + aP^Q ■ {\P - A^) ^^'^^ + ^ (fe^)' (a^ - A^)' 



= -a^ (1 - p2) (hPf (\P - (V'ip^)' + nA^ (^(") - 

[^W (A^A^,T) =n, 



(2.32) 



with initial value V'''"-' = 0, and -0^^^ — ip, as given by (2.30) 



3 properties of P^' 



Walter 



1970 



To demonstrate properties of p(") , we need a comparison principle similar to the one in 
To this end, we first state a relevant one-sided Lipschitz condition along with growth conditions. We 
require that the function g = g (A^, A^, i, v,pi,p2) satisfies the following one-sided Lipschitz condition: 
For V > w, 



g{\^ ,\P ,t,v,pi,p2) -g(A^,A^,i,w, 51,52) < c(A^,A^,i) {v - w) + di (A^,A^,t) \pi - qi\ 



+ d2 (A^A^,i) \p2-q2l 



(3.1) 



with growth conditions on c, di and ^2 given by 



(in (a^ - A^ 



In ( A^ - A^ 



<c{X^,\P,t) < K 

0<di (A^ A^, t) < X (a^ - A^) [1 + In (a^ - A') + In (a^ - A^) 
0<d2{X',XP,t) <k(xp -XP) [l + ln(A^-A') + In (a^ - A 



(3.2) 



for some constant K >0 and for all (A^, A'^,t) £ {X^ ,00) x (A^,oo) x [0,r]. 

To prove Lemma 3.2 below, as well as many of the properties of P*^"\ we rely on the following 
lemma. 

Lemma 3.1 ^C^ + <\A-B\ + ^ + P2 

Proof: It is clear that the inequality holds if A < B. For the case A> see the proof of Lemma 
4.5 inlMilevsky et all [2005]. □ 
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Lemma 3.2 Define gn, for n> 1, by 



+ aJ{l- p2) [bPf (a^ - A^) pI + n\P {v - i^i^-^))' 



(3.3) 



in which tp^" solves (2.32) with n replaced by n — 1. Then, gn satisfies the one-sided Lipschitz 



condition (3.1) on (A^jA'^jt) G (A^,cx)) x (A^,oo) x [0,r]. Furthermore, condition (3.2) holds if 



|a^''3| < K 
\aP^Q\ <K 



l + ln(A^-AM +lnfA^-A^ 



1 + In (a^ - A^) + In (a^ - A^) 



(3.4) 



for some constant K > 0. 

Proof: Suppose that v > w, then 
9n (A^, A^,t, U,Pi,P2) - 9n (A^,A^,t,w, 51,(72) 



a^'« (a^ - X') [p, q,) + a^^Q (a^ - A^) {p, - q,) nX^ [v - w) 



+ a\ {l- p2) {hP f [XP -XP] pI+ nXP {v - 



a^(l-p2)(feP)2 (^xP-X^j q^ + nXP {w - ^("-^)) 
<y''^\[x'-X')\p,-q,\+[\aP^Q\[xP-X'')+a ^1^ 6^ ( A^ - A^) ] \p2 - 92 1 

— ^nA^ — aV nXP^ (v — w) 
< |a^-«| (a^ - A') bi + [|a^'«| +av/r^6^] (a^ - A^) 1^2 -92!. 



(3.5) 



In the above series of inequahties, we use a < \/a^ < VXP and Lemma [3^" Therefore, (3.1) holds 

withc = 0, di = (a^ - A^) andds = {a^'^^l (^A^ - A^ j +a^/l-p^b 

di and ^2 satisfy condition (3.2) if (3.4 1 holds. 



P i^XP - X^j . Notice that 
□ 



Assumption 1 Henceforth, we assume that the condition (3.4) holds for rest of the paper. For later 
purpose, we also assume that a^'^ ^A^ — A^^ is Holder continuous and satisfies the following growth 
condition 

' < K 1+ flnfA^ - A^)V . (3.6) 



-P^Q IXP -XP]+aP'Q 



Theorem 3.1 Let G = {X^ ,00) x (A^,oo) x [0,T], and denote by Q the collection of functions on G 
that are twice differentiable in their first two variables and once-differentiable in their third variable. 
Define an operator C on Q by 



Cv = v, + \ {b') ' (a^ - A^) ' v,.,s +pbH^ [X' X') (a^ - A^) ^,r,. 



(3.7) 



in which gn is given by (3.3). Suppose that v,w £ Q are such that there exists a constant K >0 with 
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V < e 



{{H^'-^')rHH^^-^nr} large (ln{x' -X')y 
(in (^XP -X^^y. Then, if (a) £v > Cw on G and if {b) v {X' , X^ ,T) < w (A^A^,T) for all X^ > X' 



id A > A , then v < w on G. 



Proof: Define yi = In ^A''^ - A"^^ , j/2 = In (a^ - A^^ , and t = T - t. Write v{yi,y2,T) = 

i;(A^A^,i), etc. Therefore, v < e^(^i+^^)and w > -e^yl+y"^) for large yl + yl. Under this 
transformation, (3.7) becomes 

1 



CV^-Vr + -{¥) Vy,y, + pb^P Vy,y^ + -{bP) V y^y^ + h{yi , y2 , V , V y^ , V y^) , 



(3.8) 



in which 



Hyi,y2,T,V,Pl,P2) = - ^{b')'^Pl - ^{b^fP2 + gn{yi,y2,T,V,pi,p2), 



(3.9) 



and ti is a differentiable function defined on x [0,r]. The differential operator in (3.8) is of the 



form considered by Walter 1970 



To complete the proof, we consider the Lipschitz and growth conditions in the original variables 
A'^, A-'^, and t. From 



Walter 



1970 



we know that the conditions on h required for Walter's comparison 



principle are 

Hyi^y2,T,v,pi,p2) - Hyi^y2,r,w,qi,q2) < c{yi,y2,T){v - w) + di{yi,y2,T)\pi - qi\ 



with 



+ d2iyi,y2,T)\p2 - q2\ 

'0<c{yi,y2,T)<K{l + yf+yl), 
0<rfi(yi,y2,r) <i^(l + |?A| + |y2|), 
.0< J2(2/i,y2,r)<i^(l + |yi| + |j/2|). 



(3.10) 



(3.11) 



Under the original variables, it follows from (3.9) and (3.5) that, for v > w, 



Hyiiy2,T,v,pi,p2) - h{yi,y2,T,w,qi,q2) 
'1 



< 



,PQ 



\pi - ^1 1 



i.P.Q 



bP 



\P2 - hi- 

(3.12) 



Note that pi=e "^^pi since V-'a^ = ^ ^'i'yi'i similarly, for p2^ qi, and 52- Thus, (3.10) is satisfied with 

and d2 ^(hPf 



2 = c = 0, di = {¥f + 



3.2 



rP^Q 



a\/ \ — bP , and (3.11) is satisfied due to 



and (3.4) and the fact that b^ and bP are continuous on [0,T] and are, thus, bounded. 



Lemma 
□ 

For the remainder of this section, we apply Theorem |3.1| to investigate properties of the price 
for n pure endowment contracts. For simplicity, we will state and prove properties of t/"*^"^ and 
afterwards interpret the results in terms of P'") . 



Property 3.1 For n > 0, < 



G. 



Proof: For ease of presentation, define h by h[€) = e (~ "V^)(^ *) j^^. ^ g [0,T]. We proceed by 
induction to prove that t/'^"-' < nh on G. Note that the inequality holds for n = since V'^") = 0. For 
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n > 1, assume that ip'-^^^'^ {X' , ,t) < {n - l)h, and show that < t/A'^^X' , X^ ,t) <nh 
To apply Theorem 



3.1 



define a differential operator £ on by (|3.7|). We have Cip 



in) 



due to 



equation (2.32 1. Apply the operator C to nh to get 



C{nh) ^ ix^ - X^\ nh - {nX^ - aV^) [nh- 



< A 



n ( A^ — a\ X^ 



A^^ nh — {nX^ — a {n-{n-l))h 
(nX^ 



(3.13) 



vVnX^ 



h<0. 



Because C{nh) < = Cip'^"') and nh{T) = ip^'''\X' , X^ , T) = n, Theorem 3.1 implies that < 



ne V- ^ - ' on G. 



Similarly, we prove that ■0^"^ ^ by induction. Suppose that t/j^" > for n > 1, and show that 
■0'^") > 0. We apply the same operator C from the first part of this proof to the constant function 
on G. Because £0 (^nX^ + aVnA^^ V^""^^ > = and < n = (A^, A^, T), Theorem 



3.1 



implies that 0'-"-' > on G. 



□ 



It follows immediately from Property 3.1 that < P^"'> {r, X' , X^ ,t) < nF{r,t) for (r, A^,A^,t) G 



1 



M+ X G, in which F is the price of a T-bond with face value of $1. Thus, the price per risk — 

n 

lies between and F. This is a no-arbitrage condition since the total payoff of n pure endowments at 
time T is non-negative and is no more than $n. 

Property 3.2 For n>l, i/;'"^ > on G. 

Proof: We prove this property by induction. First, the inequality holds for n = 1 since ■0^^^ > by 



Property [3T] and = 0. For n > 2, assume that > -0^""^^ and show that 0(") > 

Define a differential operator £ on ^ by (3.7). We have that £-0^"'' = due to equation (2.32 1 



Apply the operator £ to 0^" ^\ and use the fact that 0*^" solves (2.32) with n replaced by n — 1 



(3.14) 



£^("-i) = (71 - 1)A^ (^("-1) - V--^"-')) + ^^(1 - p2) (5P)2 (a^ - A^)' (^ir'^)' 

- 0-^(1 - p2) (6P)2 (^A^ _ x'^y (^aI""'^)' + (n - 1)A^ - i^("-2))' 

> (n - 1)A^ (v/"-i) - V^("-2)^ _ a.^(n-l)A^ (i/-^"-!' - V^^"-^)^ 
= (n-l)A^-a^(n-l)A^ (^("-i) - V'""')) > 0. 

Note that the first inequality is due to the fact that \/ + B'^ < |A| + We also use the induction 



assumption that 0^" > 0^" Because 0^"-' — n > n ~ 1 — tp'^^ at i = T, Theorem 3.1 implies 
that 0(") > 0("-i) on G. □ 

We use Property |3.2| to prove Property |3.5| below; however. Property |3.2| is interesting in its own 
right because it confirms our intuition that increases with the number of policyholders. 



Property 3.3 Suppose < ai < q;2 < y A^, and let he the solution of (2.32) with a = Ui, 

for i^l,2 and for n > 0. Then, ■0(")>"i < ^(»)."2 q 
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Proof: We prove this property by induction. First, the inequahty holds for n = since 
for i = 1,2. For n > 1, assume that -i/^C"-!)'"! < ^j(«-2),a2^ ^^^j ^^^^ ^-^^^ ^{n)^ai < ^("):a2_ 

Define a differential operator £ on by (3.7) with a — ai. We have that = since 

^(n),ai golves (2.32) with a = ai. Apply the operator C to T/if")^"^ get 



(3.15) 



- ^(1 - p2) (5P)2 (a^ - A^)' +nA^ (^(«)."2 _ ^(n-l),a2)2| 



Here, we use the Lemma 



3.1 



with A = VnXP (ijA")^"'^ - , B = \^A^ (tA'"^'"^ - ^("-i),a2) 

and C = bP ^A^ — A^^ ■0^"'' "^i ^^^^ ^^'^ induction hypothesis and 012 > ai. Because -i/jt")'"! = 
^(«),ai ^ „ at t = T, Theorem El] implies that V'*")'"! < V^"^'"' on G. 



□ 



Property 3.3 shows that P*^"-' increases with the instantaneous Sharpe ratio a. The more that the 
insurance company wants to be compensated for the unhedgeable portion of the mortality risk, the 
higher it will set a. We have the following corollary of Property |3.3| 



Property 3.4 Let be the solution to ( |2.32p with a = 0. Then, for < a < \J K , V'^"^'" > 

^(").aO ^g^^ express the lower bound ^/;(")'"'^ as follows: = nip"'^ , in which ijj'^'^ 

is given by 



At — A^, Af 



A^ 



ip°'^{x\x^,t) = eQ 

and the Q-dynamics o/{A^} and {A^} follow, respectively, 

dXl = a'-'^iXl t) (Xl - X') At + b\t) (Xl X') dM//'«, 



(3.16) 



(3.17) 



dA 



f = a^^<^{Xl Af , t) (Af - A^) dt + b^it) (Af - A^) 



(3.18) 



Proof: Let ai = and a2 = « > in Property (3.3), and the inequality follows. By substituting 
a = in (2.30), the Feyman-Kac Theorem leads to the expression of -0"° in (3.16). Finally, it is 
straightforward to show that ?t,-0"° solves (2.32) with a — thus, t/j^")'"" = □ 

Note that n?/;^^-''"^ — nip'^^ is the expected number of survivors under the physical measure, so 
the lower bound of -!-P(") (as a approaches zero) is the same as the lower bound of P, namely, P-0"°. 
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Property 3.5 V^"' < on G for n > 0. 



Proof: We prove this property by induction. First, it is clear that ijj^^) = 0. For n > 1, assume 



that t/'^p < 0, and apply a modified version of Theorem 3.1 to compare i/)^"^ < and the constant 



function 0. To this end, we first differentiate ■(/'^"■''s equation, (2.32), with respect to to get an 
equation for /("^ = ^xp- 



fi 



in) 



a^'? ■ (\P - A^) + a^'Q] /(") + [a'^Q + pb' b^] (a^ - A^) + [a^'Q + (6^)'] (a^ - A^) 



2 ^'"''^ 



+ ^7, {b'r (V - y f^-l +pb'b^ A^ - y A- - A- fi'^l + ^ {bn' A- - A- fi-\^ 



(1 - p2) (5P)2 (a^ - A^) (/("))' + (1 - p2) (bPy (xP A^) /(")/, 
= —a 

In - + nXP - ^(""1)) (/(") - /("-!)) 

W (1 - p2) (6P)2 (^A-P - A^) (/("))^ + nA-P (V'C") - ^pi^-^))"^ 
/(") (A^A^,T) = 0. 



(") 



(3.19) 



Define a differential operator £ on by (|3.7p with g„ replaced by 



-[a'^Q+pb'b^] [X'-X')p, 
aP^Q + (6^)'] (a^ - A^) p2-n (^(") - t^^""!') - nA^ - /(""D) 
(1 - p2) (6^)2 (a^ - A^) «2 + (1 - p2) (6P)2 (^^P _ A^)\p2 

^(1 - p2) (6P)2 (^A-P - A^)%2 + „^P (^(«) „ ^(«-l))2 



^J{l-p^){bP f (^XP - X^y + nXP (V-H - V^(»-i))^ 



(3.20) 



From Walter 1970 , we know that we only need to verify that (3.1 1 holds for u > w = = = 52- It 



is not difficult to show that 

,/ \P ^ „ „ ^ = /\/ \P 



ff„(A^A^,^,^;,Pl,P2) - .g«(A^A^,^, 0,0,0) < 

a'^Ql+pb'b''] (X'-X') \p,\ + 



+ 1 — p^b^ V 
+ {b^f + a^/l^bP] (X^ A^) 

(3.21) 



Also, by Assumption jlj the corresponding c — a^'^ ^ 



A^-A^ 



,P,Q 



+ a^Jl- p"^ bP, di = 



ptf bP , and ^2 = |a^''5| + (t/P^ + ay^l — p^ bP in (3.21) satisfy the growth conditions in (3.2). 

Note that Cf^^'' — on G. Apply the operator C to the constant function to get £0 = 
(^nXP — aVnXP^ — — a/2y^n/ XP^ — i/;'"^^') < by the induction assumption, by 

Property Ei) and by the assumption that A^ > A^ > a"^. Since /(")(A^, A^, T) = 0, Theorem El 
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implies that /(") = V'l"^ < on G. □ 

It is intuitively pleasing that ip^j^p < because for physical survival probabilities, if the hazard rate 
increases, then the probability of surviving until time T, and thereby paying the $1 benefit, decreases. 
A related result is that P^") decreases as the risk-adjusted drift of the hazard rate, a^''^ , increases 
because the hazard rate tends to increase with its drift. 

Property 3.6 Suppose af'*^ < a^''^ on G, and let t/;^"^'"* denote the solution to (2.32 I with a^''^ — 
af^, for i = 1, 2 and for n>0. Then, ipi"-)'"-^ > o„ q 

Proof: Define a differential operator £ on by (3.7) with — af ; then, it is clear that £7/;(")^°i = 0. 
Apply this operator C to obtain 



> 0. 



(3.22) 



Since (AfA^,r) = i/-^"^'"" (AfA^,T) = n, Theo 



rem 



3.1 



implies that tjA'^')'''^ > G. 

□ 



Next, we prove the subadditivity property of P'^"'^ To that end, we use Lemma 4.10 in Milevsky 



et al. 2005 . We restate the lemma, and one can find its proof in the original paper. 



Lemma 3.3 Suppose A > C > B , B\, and C\ are constants; then, for non-negative integers m and 
n. 



(Bx + Ga) + (to + n)A^ - V^{A - C) < jBf + mB^ + JG| + nP2 + _ b). (3.23) 



Property 3.7 < f/'^"^ + V''"' for m,n>0. 

Proof: We prove this inequality by induction on to + n. When to + n = or 1, we know 
that = ip^^^ + 1//°'^ and -0(1) = V^^^^ + ■(/''■°' since V'"^ = 0. For m + n > 2, suppose that 

^e+fe) < ^j(0 -(_ ^(fe) for any non-negative integers k and I such that fc-|-/<m + ri — 1. We need to 
show that < V*"^ + V'^"^- Define ^ = ^f™) and rj = on G. The function ^ solves 

the PDE given by 



+ph^ (a^ - A^) (a^ - A^) Ca.a- + \ (bn" (a^ - A^)' Ca-.p 



-a 

v/ \P 



^{i - p^) {hPf [xp - A^)^ [i'xp)^ + (V'^"^ - 



(3.24) 



e(AfA-P,r) = m + n. 



Define a differential operator £ on C/ by (3.7) with n replaced by m -I- n. It is clear that Crj — on G. 
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Apply the operator £ to ,^ to get 
- (m + n)A^ 

+ ay(l - p2) (^A^ - A^) ' Cap + im + n)XP - 

+"-i) _ _ ^(«)^ (^jriX^ _ aV^) 



(3.25) 



{m-'rn — l) 



(m) 



nX^ 



< 0. 



To get the first inequality in (3.25), we apply Lemma 3.3 after assigning A — \/X^ — i)). 



(m) 



and Ca — 



B = (V'^'") - , c = (V'^") - v^"-'^) , Sa = bP ■ (a^ - A^) V, 

— p2 6"^ • ^A^ — A''^^ ^yp- The second inequality in (3.25) follows from the induction assumption 
^(m+n-i) ^ ^(k) ^(i) Tfj^iY^ jf^i — from the assumption that VX^ > a. Since 

£,{X',XP,T) =7] (A^, A^,r) = TO + n, Theorem O implies that < ^ on G. 



□ 



Property 3.7 states that our pricing mechanism satisfies subadditivity, p('"+") < pC™) + p("). 
This is reasonable since if subadditivity did not hold, then buyers of pure endowments could purchase 
separately and thereby save money. 

Property 3.8 Let ^(").pO be the solution to ( |232| ) with p = Oforn> 0; then, t/ji'^^-po = ^/)(")'PO(A^, t) 
is independent of X^ and solves the following PDE: 

< = -a^ibPf (a^ - A^)' (v-i"^'"")' + nA^ (^(n),pO _ ^(n-i),p0)2^ 
^^(«),po (A^,T) = 71, 

wii/i = /or n = 0. 



(3.26) 



Proof: The solution of (3.26) is independent of A-'^ and also solves (2.32) when p — 0. Uniqueness 
of the solutions of (3.26) and (2.31 ) implies that the solutions of the two PDEs are equal. □ 

When p = 0, the optimal investment in the mortality derivative is zero, as we discussed in Remark 



2.2 Also, equation (3.26) is identical to equation (4.1) of Milevsky et al. 2005 , which determines the 



price of 7i pure endowments in a market without mortality derivatives. The coincidence of the two 
results in the case of p = shows that the pricing mechanism we apply is consistent. 

It is natural to ask if the hedging will reduce the price of pure endowments. To answer this 
question, we first make an assumption on to simplify the equation for ip^"^^ as follows. 



Property 3.9 When the market price of risk for mortality is independent of X^ , then ip^"''^ — 
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■0'^")(A^,i) is also independent of \^ and solves the following PDE: 



(3.27) 



Proof: The solution of (3.271 is independent of A^ and also solves (2.32) when is independent of 



A''^. Uniqueness of the solutions of (2.321 and (3.27) implies that the solutions of the two PDEs are 
equal. □ 



Because P'"' = Ft/i^"), Property 3.9 implies that if q^' is independent of A^, then P'") is also 



independent of A^. It follows from this property and Property 3.6 that if the g'^ is independent of 
A^, then increases with increasing market price of mortality risk q^ , as one expects. 

Property 3.10 Suppose q'^ is independent of A^ and q^ < q2 ■ Let ^/jg solution of 

([3^ with q^' = q^' , for i = 1,2 and for n>0. Then, < ^(")'92' on G. 



Proof: From (l2^ we have that af'*^ > 02 , and we conclude that i/j^")^?! < -(/;("). 92 on G from 



Property |3.6| and Property |3.9 



□ 



Next, we give a condition under which hedging with mortality derivatives reduces the price of pure 
endowments. 



Theorem 3.2 Suppose q^ is independent of \^ . Letip^''^'>' denote the solution of {3.27 \ with pq^ < 



0, and let denote the solution of with p = 0. Then, " < ^/;(")'° on G 



Proof: Define a differential operator £ on by (3.7) with replaced by 



5„ {X'',t,v,p2) = ■ (a^ - A^) P2 n\P {v - ^(-D^") 
+ aJ {bPf (XP - X^YpI + nXP {v - 



(3.28) 



It is straightforward to check that the function gn in ( |3.28[ ) satisfies the one-sided Lipschitz condition 
([3T|) and the growth condtion (|3^. We have that £-0^"^'" = since solves (|3.27[) with p ^ 0. 



Apply the operator C to f/'^"-'' to get 

£,p(^^)^- =pq^'bP-(xP-x'')ti/;'J^ 



+ a^{bPf (^XP-X^y (jP'-^^J'^y + nXP (0("),- -^(n-i).o)2 
-a^{l- p2) (bPf (a^ - A^)' {'^'xp'^y + nXP (V'H-- - ^("^i).-)' 



(3.29) 



The first inequality above follows from V'l.'p' ^ 0, pq^' < 0, and Lemma 3.1 The second inequal 



ity follows by an induction step; recall that i/jC)'*^ — tp^^^^' = 0. Additionally, i/j^")' (A^,T) — 
^(n),o (A^,r), SO Theorem O implies that < V-^"^'" on G. □ 
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Remark 3.1 One can interpret the price P(") with p = as the price for which no hedging with 
the mortality derivative is allowed because the optimal investment in the mortality derivative when 
p = is 0, which follows from Property 3.8 Thus, Theorem 3.2 asserts that when pq'^ < 0, the price 



Thus, Theorem 

when hedging is allowed is less than the price with no hedging. However, if pq^ > 0, then we cannot 
conclude that hedging necessarily reduces the price of the pure endowment. We discuss this more fully 
at the end of the next section. 



4 Limiting Behavior of — p(") as n — > oo 

n 

In this section, we consider the hmiting behavior of ^P'-'^\ First, we show that the price per risk, 
kp{n) ^ decrease as n increases; that is, by increasing the number of pure endowment contracts, we 
reduce the price per contract. Then, we further explore how far ^i^^"-' decreases by determining the 
limiting value of the decreasing sequence {^^''•"''}- Surprisingly, we find in Theorem 4.1 that the 



limiting value solves a linear PDE. The proofs of most results in this section are modifications of the 

. We 



proofs given by Milevsky et al. 2005 



To prove the limiting properties of we use the Lemma 4.12 in Milevsky et al. 

restate this lemma without proof. 



2005 



Lemma 4.1 Ifn>2, and if A> C >0 and B\ are constants, then the following inequality holds 
'Bl+^C^ < V^^iA-C) + jBl + ^—[in-l)C-in-2)Af. 



77,-1 



(4.1) 



Proposition 4.1 decreases with respect to n for n > 1. 



Proof: It is sufficient to show that ^V'^"'' decreases with respect to n. Define (/)*^"^ = n'/'*'"\ £^nd 
we will show that > for n > 2 by induction. From (2.32), we deduce that 0*^"^ solves 

+pb^ b^ (a^ - (a^ - A^) ^i'j),. + \ {b^f (a^ - A^)%i'i?,P 

-A^ - (n- l)0("-i)] 

= p2) {bP f (a^ - A^) ' (0^':}) ' + ^ A^ [n<^(«) - {n - l)0("-i)] 

[^(") (A^A^,T)=1, 



(4.2) 



with 



ip, in which tp solves (2.30). 



We first show that > To this end, we define a differential operator C on Q by ( |3.7[ ) with 
gn replaced by 



52(A^ A^, t,v,pi,p2) = a^^« • (a^ - A^) pi + a^^'? • (a^ - A^) p^ - A^ {2v - ^) 



(4.3) 



It is clear that 32 satisfies conditions (3.1) and (3.2); hence, we can apply Theorem 3.1 Note that 



19 



£0(2) _ Q giijQQ 0(2) solves (4.2) with n = 2. By applying the operator C to 0(^'' — ^, we get 



(4.4) 



< = £0(2). 



Since 0(i) (A^A^,r) = 0(2) (A'f,A'f',T) 1, Theorem 3.1 implies that 0(i) > 0(2) on G. 

Assume that for n > 3, 0("-2) > 0("-i) on G, and we show that 0("-i) > 0("). We define 
differential operator £ on Q by (3.7) with g„ replaced by 



,(A^ A^, v,pi,p2) = a^-^ • (a^ - A^) pi + a^-Q ■ (a^ - A^) p^ - A^ \nv - (n - l)0("-i) 

+a^{l-p^) {bPf [\P - A^) %2 _^ ^ A-P [nv - {n - l)0(«-i)] ^ 



(4.5) 



It is clear that 5„ satisfies conditions (3.1) and (3.2); hence, we can apply Theorem 3.1 Note that 
£0(") = since 0(") solves ( |4.2[ ). Apply the operator £ to 0(" i) to get 

£0("-i) = {n- 2)A^ (0("-i) - 0("-2)) + ^^(1 - p2) (bPf (a^ - A^)' (0ir'^)' + (0("-i))' 
- aiy(l - p2) (6P)2 (^A^ - A^)^ (j^xp'^^y + ^^^^ - l)0("-i) - (n - 2)0(«-2)]2 
(n - 2)A^ - a^{n-2)\P (0("-i) - 0("-2)^ < q = £0("). 

(4.6) 



To get the first inequality in (4.6), we use Lemma 4.1 by assigning A = \/\p 0(" C ^ VXP 0(" 
and Bx = Vl-P^ bP-(^XP - A 



^ ^(" 1) We also use the induction assumption that 0(" > tp^" i). 



Additionally, 0("'i) (A^A'f',T) = 0(") (A^A'f',r) = 1, and Theorem [SjJ implies that 0("-i) > 0(") 
on G. □ 



In what follows, we answer the question inspired by Proposition 4.1 namely, what is the limit of 
the non-negative, decreasing sequence {^^'("''}? In Theorem 4.1 below, we will show the limit equals 
F/3, in which j3 = (3 {\^ ,\P ,t) denote the solution of the following PDE: 

A + a^'Q ■ (a^ - V) + [a^^Q - a^l^hP] (a^ - A^) + \ {b'Y (a^ - X')%,,,, 
+pb' bP (a^ - A^) (a^ - A^) P,,,P + I {bPf (XP A^)' /3,P,P - XPp 
-0, 

[/?(A^A^,r) =1. 

(4.7) 



By applying the Feyman-Kac Theorem to (4.7), we obtain an expression for /3 as an expectation 

/3(A^A^,i) 



Aj — A'^, Af — A^ 



(4.8) 
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in which the Q-dynamics of {A^ } and {Af } follow, respectively, 

dA,^ = a^''^{Xi, t) (\{ - A^) dt + b\t) [\l - A^) dWl' 



and 



dAf 



a^-«(A^ Af, t) - a^r^ [t]] (Af - A^) dt + (Af - A^) 



(4.9) 



(4.10) 



Here, W^, 

We begin by proving that ^P^^"^ is bounded below by F (3, and for that purpose, we need the 
following lemma. 



Lemma 4.2 The function (3 defined by (4.7) is non-increasing with respect to A 



Proof: Denote / = (3xp, and we deduce from (4.7 1 that / solves the following PDE: 



P,Q 



• A^ - A^ 



.yr^6^] (a^-a^)/,. 

, d1 „ 1 /. 2 / , r , r\ 



bP (a^ - (a^ - A^) A.,. + 1 {bPf (XP A^)' /,.,P - /3 

= 0, 

[/(A^A^,T) -0. 



(4.11) 



Define a differential operator C on by (3.7) with g„ replaced by 



5(A^ A^, t, «,pi,p2) = [a^^Q + p 6^ b^] (a^ - A^) + [a^'Q - ayi^b^] (a^ - A^) p. 



-P,Q 



A^-A^ 



v/TTp^b^-A^ 



V- 13 



(4.12) 



Because of Assumption [T] it is straightforward to check that the function g in (4.12) satisfies the 
one-sided Lipschitz condition (3.1) and the growth condtion (3.2). Because / solves (4.11), we have 
that Cf = 0. Because (3 is clearly non-negative, £0 = —(3 < 0, in which is the constant function of 
on G. Additionally, / (A^^ , X^,T) ^ 0, so Theorem [3I] implies that /S^p < 0. □ 

Lemma 4.3 For n > 1, ^-P*-"-* > F/3, in which (3 is given in (4.7) 

Proof: It is sufficient to show that ^V'^"'' > /3 on G. We prove this property by induction. First, 
for n = 1, we show that (3 < ■0'^^^ = ip. Define a differential operator £ on by (3.7) with n — 1. 
Recall that = in (lO). Since tp solves (ISlOl), Lip = 0. Also, 



£13^ a 



y (1 - p2) (feP)2 (^xP - A^) ' (31p + AP/32 _ ayr^6^ (a^ - A^) |/3^P I > = C^j. (4.13) 



Additionally, /3 (A^, A^, T) = 1 = -0 (A^ A^, T), so Theorem 3.1 implies that /3 < -0 > on G. 

For 71 > 1, assume that /3 < and show that /3 < 0^"', in which 0^") = ^^''•"^ for n > 1, as 

we defined in the proof of Proposition 4.1 Define a differential operator C by (3.7) with gn replaced 
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by gn given by (4.5). Since 0'^") solves (|4.2|), — 0. By applying this operator on /3, we get 



- av/r^ fo^ (a^ - A^) I + A^ [{n - - (n - 

>0 = 



(4.14) 



Also, P{\^,XP,T) = (A^A^,r) = 1; thus, Theorem 3.1 implies that P < = ^V^"^ on G. 
□ 

Next, we show that hm„_^oo -P^"'^ = Pp. To this end, we need some auxiliary results. First, we 
prove that V'^"^ is bounded from above by 7^"^ = 7'^"^(A^, A^, t) for n > 0, in which the function 7*^"^ 
solves the following PDE: 



'7(") + a^^Q ■ (a^ - A^) 7i"' + - a^T^b^] (a^ - A^) + \ {b^f (a^ - A^)' 7^'; 



+pb^ bP (a^ - (a^ - A^) 7i"i. + ^ {bn 

- (nA^ - aV^) (7(") - ^ 0, 

[7("' (A^A^,r)=n, 



in which 7*^°' = 



A^-A^ 



■ (n) 



(4.15) 



Lemma 4.4 The function ^^'^^ given by (4.15) is non-increasing with respect to \P , and^^"^"^ > 7*^" 
for n > 1 on G. 



Proof: The proof that 7^"^ < is similar to the proof that i/j^'p < in Property 3.5 Also, the proof 
that 7^"-' > 7("^i) is similar to the proof that ?/'*•"'' > i/'''""^'' in Property 3.2 Therefore, we omit the 
details of the proof. □ 

Lemma 4.5 For n > 0, 7^") > V^"^ on G. 

Proof: We prove this lemma by induction. For n — 0, we have 7*^"^ = V'^"'' — 0. Assume that for 
n > 1, we have 7("-i) > ■^/'^""^^ and show that 7^") > V'"^- For this purpose, define a differential 
operator £ on t/ by (3.7). Then, Ctp — 0, and 

£7("' = av/r^6^ (a^ - A^) 7^"; + (nA^ - aV^) (j'^"^ - 7("-i)) - nX^ (j'^"^ - ^^^-^^') 



in) 



< 



(4.16) 



The first inequality above is due to the fact that 7^"^ < 0, that 7^") > 7("-i) > and that 

+ 52 < \A\ + \B\. Additionally, we have that 7(") (A^,A^,T) = V^"' (A^AP,r) = n; then. 
Theorem [3T]imphes that 7^") > on G. □ 

Next, we prove the main result of this section. 
1 



Theorem 4.1 lim„ 



P(")(r, A^A-P,t) = F{r,t)p{X',XP,t) on G. 
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Proof: By Lemmas 4.3 and 4.5 it is sufficient to show that hm„_j.oo (^7^"^ ~ = since ^7^"^— /3 > 
_ ^ > 0. For n > 1, define T^") on G by T^") = i7(") - /3, so we just need to prove that 
fim„^oor(") = 0. For n>l, the function F^") solves the following PDF: 



- ) F^'i.' 



+ -Abr(>^'->^') A^ - AM A- - An r[-\^ + t (b^y _ A 



.(") 



^-(n-1) A^-a 



p(n-l) 



F(") (A^A^,T) =0, 

(4.17) 

with < F*^^) = 7^^) — /? < 1 on G. By applying the Feyman-Kac Theorem to (4.171, we obtain the 
following expression for F(") in terms of f("-i): 



F(")(A^A^,i) = aE^ 
+ {n- 1)E' 



du 



ds 



Ar-a 



Aj — A^, Af — A^ 



(4.18) 



in which the Q-dynamics of {A^ } and {Af } follow, respectively, equations (4.91 and (4.10). 

Suppose F^""^) < Kn-i on G for some n > 2 and for some constant Kn-2 > 0. Note that /3 < 1 
on G, so we get the following inequality: 



F(")(A^A^,^) <a¥}^ 
+ {n- l)i^„_i ¥P 



Af -//(riA^;-a^A^)du 



ds 



Af — a\ 



n 



du 



ds 



(4.19) 



Equivalently, we can write the inequality (4.191 as 



r("HAfA^,t) < (Af A^,t) + I^/f„_ii?(") (Af A^,i) 



in which the functions and B*^"-* are defined as 

i-T 



(Af A^,t) =aEQ 

and 

(Af A^,i) =£"5 



ds 



At — A^, Af — A^ 



(4.20) 



(4.21) 



?iAf -aynAf^ e ■'''°( 



ds 



A* — A'^, Af — A'^ 



(4.22) 



After the next two lemmas that give us bounds on A*^"^ and i3("f respectively, we finish the proof of 
Theorem O □ 
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Lemma 4.6 Forn>2, < J = 



V2 



on G, in which A*^"^ is defined in (4.211 



Proof: By the Feyman-Kac Theorem, A^") in (|4.21|) solves the following PDE 



/((") 



' + a^'Q • (a^ - A^) ^-J) + [a^-'? - a^/l^bP] (a^ - A^) A^^'J + ^ {h'f (a^ - A^) 
+pb' hP (\' - A^) (a^ - A^) A^;^}^, + \ {bPy (a^ - A^)' - (nXP - aV^) 

= ~an\/ A^, 

(A^A■P,^) =0. 

(4.23) 



For n > 2, we define a differential operator £ by (3.7 1 with g„ replaced by 



„(A^A^,^,^;,Pl,P2) = a'''^• fA^-A^)pi+ [a^'«-av/r^6^1 f A^ - A^) p2 



Since ^„ satisfies conditions (3.1) and (3.2 1, we can apply Theorem 



3.1 



(4.24) 



It is clear that LA^'^'i = 0, 



and by applying the operator £ to J, the function that is identically equal to J, we get 
£J = - (n\P - a%/rJ^) J + an\fW < = £A("). 



Since A^") (A^,A^,T) = < J, Theorem 3.1 implies that < J on G 



(4.25) 

□ 



Lemma 4.7 For n>2, S^") < 1 on G, m which B^"'> is defined in (|4.22|) 



Proof: By the Feyman-Kac Theorem, S'") in (4.22 1 solves the following PDE 

+ a^'Q ■ (a^ - A^) B^;^ + [aP^Q - av/T^fe^] (a^ - A^) B^;'J + ^ {h'f (a^ - A^) 



B 



(n) 



-pb' bP (a^ - A^) (a^ - A^) B^'J^^p + ^ {bPf (\P - A^)' B^';\p - (n\P - aVnA^) 



= — yiXP — a\/n\P 

(A^A-p,r) = 0. 



For n > 2, we define a differential operator £ on by (3.7) with replaced by 



g„(A^A^,i,t,,pi,P2) = a^''^ • - A') PI + -«yr^&^] - A^) 



(4.26) 



(4.27) 



Since 5„ satisfies conditions (3.1 ) and (3.2 ), we can apply Theorem 3.1 It is clear that £_B*^"^ = 0, and 
by applying the operator C to 1, we get £1 = — {nXP — a^/ nX^^ + ^nA''^ — a\/ nX^^ = = C B . 
Since (A^, A^, T) = < 1, Theorem ET] implies that < 1 on G. □ 



End of Proof of Theorem \4-l\ By Lemmas |4.6| and |4.7[ we get the following result: for n > 2, if 

r("-i) < Kn-i, then 



r(") <K ^ — + "~ V 1 
nyi'^ n, 



(4.28) 
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with Ki = 1. Define L„ = nKn and note that i„ = in-i 



J 



for n > 2. It follows that 



J 

«=2 



dx 



<l + 2JVf^, n>2, 



which implies that on G, 



r(") < i4:„ < - + n > 1. 



(4.29) 



(4.30) 



lim„_j.oo ^ + "7= — 0; therefore, F^") converges to uniformly on G as n goes to infinity. In other 



words, lim„^oo ^P^"'> = F /3 on G. 

We end this section with some properties of /3 with the goal of determining the effect of p on /3. 

Property 4.1 If q'^ is independent of \^ , then (3 ~ /3(A^,i) is independent of and solves the 
following PDE: 



0, 



(4.31) 



/3(A^,T) =1. 

Proof: The solution of (4.31) is independent of A^ and also solves (4.7 1 when q'^' is independent of 



A''^. Uniqueness of the solutions of (4.7) and (4.31 ) implies that solutions of the two PDEs are equal. 

□ 



Theorem 4.2 Suppose q^ is independent of , and define a ^ — pq^ + a-^/l — 



Let 



denote the solution of (4.31) with a — hi, fori — 1,2. Then, (5°'^ > fS""^ on G if ai < 0,2- 



Proof: Define a differential operator £ on by (3.7) with g„ replaced by 



g(A^z;,p)=ai- (a^-A^)p-A^z;. 



(4.32) 



It is straightforward to check that the function g in (4.32) satisfies the one-sided Lipschitz condition 



(3.1) and the growth condition (3.2). Since I3°-^ solves (4.31) with a ~ ai, we have that Cf3°-^ = 



Apply this operator on to obtain 



£(3'^' = (ai - aa) (a^ - A^) /^^^ > = CP^'K 



Since (3^^ (X^,T) = (A^,T) = 1, Theorem 3.1 implies that {3^^ > l3^^ on G 



(4.33) 

□ 



Remark 4.1 When p ~ 1, namely the the insured individuals and the reference population face 
the same uncertainty in their respective hazard rates, the limiting price per contract is reduced by 
hedging when q^ is less than the pre-specified instantaneous Sharpe ratio a. Indeed, the drift a = 



a 



from Theorem 



4-2 equals a^ — q^ b^ when p = 1- Also, the effect of not 
allowing hedging can be achieved by setting p ^ throughout our work, as discussed in Remark \2.2\ 
in that case, the drift a becomes a^ — ab^. Thus, according to Theorem 4-2 the limiting price per 



contract is reduced when hedging is allowed if q^ < a. 

In other words, hedging with mortality derivative benefits the insured, through a reduced price, 
when the market price of mortality risk is lower than that required by the insurance company. In this 
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limiting case, the risks inherent in the contract can be fully hedged using the interest rate derivative 
and the mortality derivative. Indeed, the variance of the hedging portfolio goes to as n goes to infinity 
when p = I- Refer to Remark \2.1\ in which we discuss the mortality risk in the single-life case. So, as 
n — > oo, the risk coming from the timing of the deaths disappears; compare with (2.27). 



The price of the contract is reduced by transferring the mortality risk to a counterparty who requires 
a lower compensation for the risk than the insurance company does. By contrast, for a single pure 
endowment contract, the volatility in the contract due to the uncertainty of the individual's time of 
death is not hedgeable with mortality derivatives even when p = 1. In the single-life case, even if 
is less than a, hedging does not guarantee a reduction of the contract price. 



Corollary 4.1 Suppose q^ is independent of \^ , and let denote the solution of (4.311 with 



qP = at, for i = 1,2. Then, > 13"-'- on G if ai < 02. 

The result above is consistent with our intuition. Indeed, with a higher drift on the hazard rate, 
the individual is less likely to survive to time T, and, consequently, the (limiting) value of the pure 
endowment contract is lower. 



Corollary 4.2 Suppose q^ is independent of , and let f3f' denote the solution of ( 4.31[ ) with p — Pi 
fori = 1,2. Then, < (3p^ on G if pi q^' + - p\ < p2 q^' + ay^l - p\ for all t G [0, T] . 



Remark 4.2 A natural question that follows from Corollary is when is f{p,t) = pq (t) + 



ay^l — decreasing with respect to p for t G [0,r]? Suppose that p > 0, which is what one ex- 
pects between the insured and reference populations. If f is decreasing with respect to p > 0, then 
greater positive correlation will lead to a lower per- contract price, an intuitively pleasing result. It is 
straightforward to show that f decreases with respect to p if and only if 



P > 



(4.34) 



This inequality holds automatically if q^ < 0, that is, if the mortality derivative is a so-called natural 
hedge, which we discuss more fully in Remark 4.3 below. When q^ > 0, it holds for p in a neighborhood 
ofl. 

We have the following special case of Corollary |4.2[ 

Corollary 4.3 Suppose q^ is independent of \^ . If pq'^ + a-\/l — p^ < a, then the limiting price 
per risk in which hedging is allowed is less than the limiting price with no hedging (p = 0). 

Remark 4.3 In particular, when q^ is negative (and p is positive), the unit price of the contract 
is reduced by hedging, as demonstrated in Corollary \4.3\ Since the correlation is usually positive, a 
mortality derivative with a negative market price of risk q^ , that is, a natural hedge, is preferred. An 
example of a natural hedge is life insurance, as discussed in Young f200^ , although strictly speaking 
this insurance product is not a mortality derivative traded in the financial market. Both \Bayraktar\ 
and Young \2007a ' and Cox and Lin 12007^ proposed hedging pure endowment or life annuity contracts 



with life insurance. 
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5 Numerical Example 



In this section, we demonstrate our result with numerical examples. We assume that the risk-free rate 
of return r is constant and focus on the effect of the correlation p and the market price of mortality risk 
. We also assume that the market price of mortality risk is constant, and, thereby, is automatically 
independent of A^. In this case, P*^"-' and lim„^oo do not depend on A'^, as we prove in Properties 



3.9 and 4.1 Moreover, we assume that the hazard rate A^ follows the process in \2.2\ with and 



constant. We compute the price for a single contract, P(r, A-'^, t) — e~^^"^~*^ 1 1)^ a-nd the limiting 
price per contract for arbitrarily many insureds, hm„_>.oo ^P'^^\r,\^ ,t) = e^''^-^^*) /3(A^, t), and we 
use the following parameter values: 

• The pure endowment contract matures in T = 10 years. 

• The constant riskless rate of return is r = 0.04. 

• The drift of the hazard rate is — 0.04. 

• The volatility of the hazard rate is — 0.1. 

• The minimum hazard rate of the insured individuals is A^ — 0.02. 

• The risk parameter is a — 0.1. 



See Section 
In Figure 



mfo: 
■e[zl 



for the algorithm that we use to compute ip and j3. 

for a variety of values of the market price of mortality risk , we present the price 
of a single-life contract P and the limiting price per contract lim„_j.oo ^-P''"'- It follows from Theorem 
|4.2| that, given a positive correlation p > 0, the limiting unit price of a pure endowment is greater with 
a greater market price of mortality risk q^ , and the second set of graphs in Figure [t] demonstrates 
this result. Notice that the price of the unhedged contract is the price with p = 0. Since in the pricing 
mechanism, we hedge the volatility with the mortality derivative as much as possible to reduce the 
variability of our hedging portfolio, a large value of could lead to a higher contract price than that 
of an unhedged one. Observe this numerically in graphs in Figure [Tj 

we concluded that \i q^ < a and if p = 1, then the limiting price per contract is 



In Remark 



4.1 



less than the limiting price per contract of an unhedged portfolio of pure endowments. This result is 
supported by our numerical work; indeed, the curve for q^ = 0.15 in the second set of graphs lies 
above the unhedged price of approximately 0.343, the price when p = 0. 

In that same remark, we noted that for a single-life contract, we cannot conclude that the price 
with hedging will be smaller than the price without hedging, even when q^ < a. This conclusion is 
also supported by our numerical work; indeed, the curve for q^' = 0.09 in the first set of graphs lies 
above the unhedged price of approximately 0.435, the price when p — 0. 

Figure [7] also demonstrates the relation between the unit price of a contract and the correlation 
p. Take the limiting price per contract lini„_j.oo '^P^^^ with q^ 0.05, for example. When p = 1, 
hedging is preferred to not hedging, in terms of reducing the price of the contract. By contrast, 
when p < 1, that is, the two mortality rates A^ and A^ are not perfectly correlated, hedging may 
increase the unit price of the contract such as the case when p = 0.8. This observation indicates 
that the population basis risk, which is the risk due to the mismatch of the insured population and 
the reference population, diminishes the effectiveness of hedging. This mismatch, or equivalently, a 



correlation p < 1, may lead to a higher unit price for the hedged contract. See Coughlan et al. 2007 
for discussion of population basis risk. 
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6 Conclusion 



In this paper, we developed a pricing mechanism for pure endowments, assuming that the issuing 
company hedges its pure endowment risk with bonds and mortahty derivatives, and requires com- 
pensation for the unhedged part of mortality risk in the form of a pre-specified instantaneous Sharpe 
ratio. In our model, we took the hazard rates of the insured population and reference population, as 
well as the interest rates, to be stochastic. We derived the pricing formulae for the hedged contracts 
on single life and on multiple conditionally independent lives. Also, we obtained the pricing formula 
for the limiting price per pure endowment contract as the number of the insureds in the portfolio goes 
to infinity. In each case, the price solves a PDE, and we analyzed these PDE and thereby determined 
properties of the prices of the hedged pure endowments. The limiting price per contract solves a linear 
PDE and represent this value as an expectation with respect to an equivalent martingale measure. 
We noted that, in the limiting case, the mortality risk inherent in the pure endowment is fully hedged 
by the mortality derivative when the correlation between the two hazard rates and is 1. 

To investigate the factors that affect the effectiveness of hedging, we devoted our attention to the 
market price of the reference mortality risk and the correlation p between A^ and A^ . Since the 
correlation p is more likely to be positive in reality, we focused on the case for which p > during our 
discussion (and especially in our numerical work) and assumed that the market price of the mortality 
risk is independent of A^. We found that hedging with a mortality derivative requiring a negative 
market price of mortality risk always reduces the price of the contract. This result is consistent with 
the conclusions in Bayraktar and Young 2007a and Cox and Lin 2007 that hedging pure endowments 
(or life annuities) with life insurance reduces the price of the former. 

For the limiting case, we reached a more straightforward conclusion, as we discuss in Remark 
Specifically, if p = 1, the condition that < a guarantees a reduction in the per-contract price 



4.1 



through hedging. However, ii p < 1, it is possible that hedging with the mortality derivatives increases 
the price of the contract even if this condition is satisfied. This result reflects the significance of p on 
the effectiveness of hedging. We also found that, in our numerical work, hedging with the mortality 
derivatives is less effective in reducing the variance of the hedging portfolio for pure endowments of a 
finite number of individuals. 

Our results suggest that, to make it efficient for underwriters to hedge mortality risk and thereby 
benefit the insured, transparent design of mortality indices and mortality derivatives is essential. 
Reducing the market price of the mortality risk q^ is also critical. Therefore, it is important to build 
up a liquid mortality market and provide more flexible mortality- linked securities in order to reduce 

In our paper, we only investigated the prices of pure endowments and assumed that the mortality 
derivative is a g-forward. However, we believe that the main qualitative insights will hold in general. 



7 Appendix 

In this section, we present an algorithm for numerically computing ip. Recall that in our numerical 
example, we assume that q^ , the market price of mortality risk A^, is a constant, as well as and 
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. Then, equation (2.30) becomes 



+ [a^ - pq^'b^] (a^ - A^) + I {b^f (a^ - A^)%,P,P - A^V 

= -a^{l-p^) [bPf (\P - A^) ' + A^V', ^^-^^ 
[V'(A^,r)=l. 

Next, we describe our numerical scheme to compute if). 
Transformation Define t = T — t, y = In ( A^ — A^ J , and ip (y, t) = ij) (A-^, . By (7.1 ), '4> solves 



(7.2) 



^^(y,o) = 1, 

in which a = — pq^ ^ \ i^^)^ ■ 

Boundary Condtions While equation ( |7.2[ ) for ip is defined in the domain K x [0,T], we solve it 
numerically in the domain [—M,M] x [0,r] such that e~^^ is approximately zero. Therefore, 
we require boundary conditions at y = ±M . 



1. If Af A^, then Af A^ for all s e [t, T]. From equation ( [7\T| ), we have that (a^, = 
exp < — I A^ — ay ] {T — t) >. Thus, it is reasonable to set the boundary condition at 



y = -M to be V {-M, t) = exp |- - a \l k j t 

2. If A^ is very large, we expect the individual to die immediately, so the value of the pure 
endowment is approximately 0. Thus, we set the boundary condition at y = M to be 
ip {M, r) = 0. 



Finite Difference Scheme We discretize the differential equation ( 7.2 ) and get a corresponding 
difference equation as follows: 

1. Choose the step sizes of y and t as h and k, respectively, so that / — 2M/h and J — T/k 
are integers. 

2. Define yi = —M + ih, Tj = jk, and tpij — ijj (yi, Tj), for i = 0, 1, . . . , / and j = 0, 1, . . . , J. 

3. We use a backward difference in time, central differences in space, and a forward difference 
for the square-root term. Therefore, we have the following expressions: 



'4'y iVi^Tj) 



■pi+l.j + l - ■Pt-l,]+l ^^(,2\ 



(7.3) 
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Also, for the non- linear term in (7.2 1, we have 



\ 



(1 - P-) ibn' I ^'^^-^^^-"^ ) + (e- + + Oih). 



(7.4) 



Therefore, we approximate (7.2) to order 0{k + h) with the following difference equation: 

(7.5) 



- jii.] _ . V't+ij+i ~ ,\(iP\'^ ipi+i,j+i - '2'4'i.j+i + V^i-ij+i 



in which 



(l-p2)(foP)2 



"0^+1 J - 



2/i 



(7.6) 



If we define a^^ - {b^Y 2^2'^^!+ (^'^)^ 7? + ^ A^' '^^^d -a^ - {b^Y then 



(7.5 1 becomes 



aV^i-ij+i + (6 + fce^') TAi^j+i + c-0i+ij+i = V^ij + afc^ij, (7.7) 
for i = 1, 2, ...,/— 1 and j = 0, 1, . . . , J — 1, with the following boundary conditions: 



[b + /ce^i) V'lj+i + c'02,j+i = V'lj + akAij — ae 
a'(/>/_2j+i + {b + key'-^)ipi-i,j+i = '^/-ij +afc^/-i,j. 



(7.: 



for J = 0, 1, . . . , J — 1. It is convenient to write equations (7.7)-(7.8l in matrix form as 

(7.9) 



M*j+i = *j + akAj - 



,0,...,0 



for j ~ 0, 1, . . . , J — 1, in which the superscript t represents matrix transpose. In the 



equation above, '4'j 



V'lj,'02,i, ■ • ■,i>I-\,3 



and Aj = [^i.jj ^2j": ■ • ■ j ^/-ij]* with A 



defined in (7.6). The matrix M is a tri-diagonal matrix with the sub-diagonal identically 



a, with the main diagonal 6 + fce^^ , b + fce^^ , . . . , 6 + ke^'^^ , and with the super-diagonal 
identically c. 

4. Begin with the initial condition ipi^o = 1, for i = 1, 2, ...,/— 1, or equivalently, 'S'o — 1, 



in which 1 is an (J — 1) X 1 column vector of Is. Then, solve (7.9) repeatedly for j = 
0, 1, . . . , J — 1 until we reach ^j. 

One can modify this algorithm to compute V'*-"'' for any n > 1 and the limiting result /3 = lim„_i.oo V''"''- 
Computing the latter is particularly straightforward because /3 solves a linear PDE. 
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